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Abstract 

We study orthogonal matrices whose elements have moduli < 1. This paper shows that 
the existence of two such families of matrices is equivalent. Specifically we show that the 
existence of an Hadamard matrix of order 4 1 is equivalent to the existence of a 2-level 
Cretan-Mersenne matrix of order 4t- 1. 
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1. Introduction 

An application in image processing led to the search for orthogonal matrices, all of 
whose elements have modulus < 1 and which have maximal or high determinant. 

Cretan matrices were first discussed during a conference in Crete in 2014 by N. 
A. Balonin, M. B. Sergeev and colleagues of the Saint Petersburg State University of 
Aerospace Instrumentation but were well known using Heritage names , [1, 2, 3, 4]. This 
paper follows closely the joint work of N. A. Balonin, Jennifer Seberry and M. B. Sergeev 
[5, 6, 7], 

In this and future papers we use some names, definitions, notations differently to how 
they have been have in the past [2]. This we hope, will cause less confusion, bring our 
nomenclature closer to common usage, conform for mathematical purists and clarify the 
similarities and differences between some matrices. We have chosen to use the word level, 
instead of value for the entries of a Cretan matrix, to conform to earlier writings [2, 3, 4]. 

We know of no references where the Co-Existence of Hadamard Matrices and 
Cretan(4f - 1,2)-Mersenne Matrices Theorem is stated. 

1.1. Preliminary Definitions 

The absolute value of the determinant of any matrix is not altered by 1) interchanging 
any two rows, 2) interchanging any two columns, and/or 3) multiplying any row/or column 
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by -1. These equivalence operations are called Hadamard equivalence operations. So the 
absolute value of the determinant of any matrix is not altered by the use of Hadamard 
equivalence operation. 

Although it is not the definition used by purists we use orthogonal matrix as below. 
Write I n for the identity matrix of order n and let w be a constant. When a matrix S is 
written in the following form 

x a ... a 
a 



a 


B is said to be the core of S and the a’s are the borders of B in S. The variable x is set 
as zero for conference matrices, cr for Hadamard matrices and a for Fermat matrices. 

Definition 1 (Orthogonal Matrix, Hadamard Matrices, Cretan Matrix). An 

orthogonal matrix , S = ( Sij ) of order n, is square and has the modulus of all its entries 
< 1, and satisfies SS T = uil n for n = 1,2,4 1. 

An Hadamard matrix of order has entries ±1 and satisfies HH T = nl n for n = 1, 2, 4 1, 
t > 0 an integer. Any Hadamard matrix can be put into normalized form , that is having 
the first row and column all plus Is using Hadamard equivalence operations: that is it 
can be written with a core. 

A Cretan matrix , S, of order v has entries with modulus < 1 and at least one element 
in each row and column must be 1, and which satisfies SS T = ujI v . A Cretan(n]T\uj ) 
matrix, or CM(iv,t\uS) has r levels or values for its entries. 

In this work we will only use orthogonal to refer to matrices comprising real elements 
with modulus < 1, where at least one entry in each row and column must be one. 

Hadamard matrices and weighing matrices are the best known of these orthogonal 
matrices. We refer to [6, 8, 9, 10] for more definitions. We recall Barba [11] showed that 
for matrices whose entries have modulus < 1, B , of order n 

det B < \/2n - l(n - l)“s~ or asymptotically s» 0.858(n) ¥ . (1) 

Wotjas [12] showed that for matrices whose entries have modulus < 1, B , of order 
n = 2 (mod 4) we have 

det B < 2(n - l)(n - 2) or asymptotically 0.736(n) ^ . 

More details of Cretan matrices can be found in [13]. 

When Hadamard introduced his famous inequality [9], for matrices with moduli < 1 it 
was noticed that Hadamard matrices, which are orthogonal and with entries ±1 satisfied 
the equality of Hadamard’s inequality. 

However the inequality applies to other matrices and orders with entries on the 
unit disk, these have been named “Fermat”, Mersenne and Euler matrices for different 
congruence classes (mod 4), [6, 13]. 

We emphasize: in general a Cretan(n, r) or CM(n,r) or S, has r levels, they are 
made from r-variable orthogonal matrices by replacing the variables by appropriate real 
numbers with moduli < 1, where at least one entry in each row and column is 1. After 
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this stage, r-variable orthogonal matrices and Cretan(n,r), are used, loosely to denote 
one-the-other. 

For 2-level Cretan matrices we will denote the levels/values by x, y where 0 < |y| < x = 1. 
We also use the notations Cretan(v), Cretan(v)-SBIBD and Cretan-SBIBD for Cretan 
matrices of 2-levels and order v constructed using SBIBDs. 

Definition 2 (SBIBD and Incidence Matrix). For the purposes of this paper we 
will consider an SBIBD(v, k, A), B, to be a v x v matrix, with entries 0 and 1, k ones per 
row and column, and the inner product of distinct pairs of rows and/or columns to be A. 
This is called the incidence matrix of the SBIBD. For these matrices X(v - 1) = k(k - 1). 

For every SBIBD(v,k, A) there is a complementary SBIBD(v,v - k,v -2k + A). One 
can be made from the other by interchanging the 0’s of one with the l’s of the other. The 
usual use SBIBD convention that v > 2k and k > 2A is followed. 

A combinatorial trick allows us to say that any matrix of order n satisfying AAA = 
a I + bJ has determinant ( \/a + nb ) a~. 

Writing A for the incidence matrix of the SBIBD{v , k, A), which has entries 0 and 1, 
and B = 2A - J, J the matrix of all ones, for the ±1 it forms we have: 


AA T = (k - A) I + A J AJ = kJ and 

BB J = 2(fc - A)/ + (v - 2(k - A)) J BJ=(2k-v)J. (2) 


Thus 

det(A) = k(k - A)~ and det(B) = \Jk 2 + (v - k) 2 )(2(v - k)g~ . 

We now define our important concepts the orthogonality equation , the radius equa- 
tion(s), the characteristic equation(s) and the weight of our matrices. 

Definition 3 (Orthogonality equation, radius equation(s), characteristic equa¬ 
tion^), weight). Consider the matrix S = (s,j) comprising the variables x and y. 

The matrix orthogonality equation 

S T S = SS T = wl n 


yields two types of equations: the n equations which arise from taking the inner product 
of each row/column with itself (which leads to the diagonal elements of u>I n being ui) are 
called radius equation(s), g(x,y) = co. and the n 2 - n equations, f(x,y) = 0, which arise 
from taking inner products of distinct rows of S (which leads to the zero off diagonal 
elements of uil n ) are called characteristic equation(s). The orthogonality equation is 
T.]=i s 2 j = u). oj is called the weight of S. □ 

Example 1. We consider the 2-variable S matrix given by 


x 


S = 


V 

y 

y 

.y 


y y y y 

x y y y 
y x y y 
y y x y 
y y y x_ 
3 



By definition, in order to become an orthogonal matrix, it must satisfy the orthogonality 
equation, SS T = uil, the radius and the characteristic equations, so we have 

x 2 + 4 y 2 = w, 2 xy + 3 y 2 = 0 . 

To make a Cretan(5;2;-y) we force x = 1, (since we require that at least one entry 
per row/column is 1), and the characteristic equation gives y = -|. Hence oj = 3|. The 
determinant is (yy) = = 20.286. Thus we now have anS= Cretan{ 5; 2; 20.286). □ 

1.2. Notation Transitions 

In transiting from one mother tongue to another (Russian to English and English to 
Russian) and from previous to newer usage, some words reoccur: we need a shorthand. 
To simplify references we note: 


Heritage Usage 

Cretan Matrix 

References 

Fermat 

Cretan(4t+1) 

[2, 6, 14, 3] 

Hadamard 

Cretan(4t) 

[5, 4, 6, 9, 8] 

Mersenne 

Cretan(4t-1) 

[1, 2, 4, 6, 15, 16] 

Euler 

Cretan(4t-2) 

[1, 3, 4], 


Table 1: Cretan and Heritage Names 


Usage 

Dual Usage 

Heritage name 

CM(4t+l) 

CM(4t) 

CM(4t-l) 

CM(4t-2) 

CM(4t-3) 

Core of CM {At + 2) 
Core of CM {At +1) 
Core of CM {At) 

Core of CM {At - 1) 
Core of CM {At - 2) 

“Fermat” matrices 
Hadamard matrix 
Mersenne matrix 
Euler matrix 
“Fermat” matrices 


Table 2: Cretan and Core Names 


2. Preliminary Results and Hadamard Mathematical Foundations 

As an historical note we point out that J. A. Todd’s [17] article showing the relationship 
of SBIBD{At - 1,2f - l,t - 1) and Hadamard matrices of order At appeared in the same 
issue of the Journal of Mathematics and Physics as the famous paper by R. E. A. C. 
Paley [18] using Legendre symbols to construct orthogonal matrices. 

Example 2. An SBIBD(7, 4,2) = B and its complementary SBIBD{ 7,3,1) can be writ¬ 
ten with incidence matrices: they are still complementary if permutations of rows/columns 
are applied to one and other permutations of rows/columns to the other. This is because 
if P and Q are permutation matrices PBQ is equivalent to the SBIBD(7,A,2). 
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'1 1 1 0 1 0 o' 

0 1110 10 
0 0 1110 1 
SBIBD{ 7,4,2)= 10 0 1110 
0 10 0 111 
10 10 0 11 
110 10 0 1 

'0 1 1 0 1 0 o' 

0 0 110 10 
0 0 0 1 1 0 1 

SBIBD{ 7,3,1)= 1 0 0 0 1 1 0 

0 1 0 0 0 1 1 

1 0 1 0 0 0 1 

_1 1 0 1 0 0 0_ 

The second incidence matrix is still a complement of the incidence matrix of the first 
SBIBD even after permutations of its rows and/or columns have been performed. 

To prepare to make Cretan matrices from SBIBD(4:t - 1,2 1— 1 ,t- 1) we first put 
the two SBIBD in variable forms: or SBIBD(7, 4,2) circ (x,x,x,y,x,y,y) and for the 
complementary SBIBD(7, 3,1), circ (y,x,x,y,x,y,y), for 

x x x y x y y 

y x x x y x y 

y y x x x y x 

SBIBD{ 7,4,2)= x y y x x x y 

y x y y x x x 

x y x y y x x 

x x y x y y x 

■y x x y x y y 

y y x x y x y 

y y y x x y x 

SBIBD{ 7,3,1)= x y y y x x y 

y x y y y x x 

x y x y y y x 

x x y x y y y 

We could just have exchanged x and y in the definition circ(a;, x, x, y, x, y, y) to get 
the second SBIBD , but we wanted to illustrate that there are many other possibilities 
for the first row of the second SBIBD. 

Consider the 2-variable SBIBD{7 ,4,2): it has characteristic equation 2a; 2 +4a ry+y 2 = 0, 
and radius equation u = 4a: 2 + 3 y 2 . det(S') = w 5 . This is an CM( 7,2). To make a Cretan 
matrix we now set x = 1 and solve the characteristic equation to find y in terms of x. This 
value/level is then used to give the Cretan(7; 2; 5.0294). Thus the principal solution has 

x = 1, y = -2 + \/2, ui = 4a; 2 + 3 y 2 = 5.0294 . 
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(a) the principal solution (b) the complementary solution 


Figure 1: Orthogonal matrices for order 7: Balonin-Sergeev Family 


The SBIBD(7, 3,1) with characteristic equation x 2 +Axy+2y 2 = 0, and radius equation 
to = 3a; 2 + 4 y 2 , det(S') = o7 , (smaller than above values replacing y ), has a feasible solution 
for the characteristic equation. Proceeding as before we obtain a second solution 

-2 + \/2 o n 

x = 1, y = ---, w = 3x 2 + 4 y - 3.3431 . 

Thus it gives a Cretan(7; 2; 3.3431) matrix. 

The two determinants are 285.31 and 69.319 respectively. 

Loosely we write Cretan{7\2 : 5.2904) - SBIBD, (2-variable orthogonal SBIBD 
(7,4,2)), and obtain the following 2-level Cretan-Mersenne matrix: 


1 

1 

1 

-2 + 72 

1 

-2 + 72 

-2 + 77 

-2 + 72 

1 

1 

1 

-2 + 72 

1 

-2 + 72 

-2 + 72 

-2+72 

1 

1 

1 

-2+72 

1 

1 

-2 + 72 

-2 + 72 

1 

1 

1 

-2 + 72 

-2 + 72 

1 

-2 + 72 

-2+72 

1 

1 

1 

1 

-2 + 72 

1 

-2 + 72 

-2 + 72 

1 

1 

1 

1 

-2 + 72 

1 

-2 + 72 

-2 + 72 

1 


3. Cretan-SBIBD(v,2) Theorem: The Cretan 2-level Matrices from SBIBD 
Theorem: 

Although our main theorem appeared as a corollary in [13] it is in fact worthy of being 
a theorem in its own right. The paper [13] gives many relevant details. We write 
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Theorem 1 (The Cretan-SBIBD(v,2) Theorem: The Cretan 2-level Matrices 
from SBIBD Theorem). Whenever there exists an SBIBD(v,k,X ) there exists a 
2-level Cretan-SBIBD(v,2), or S or CM, as follows, 

• Cretan(y; 2; kx 2 + (v - k)y 2 \ y,x;determinant) ; 

• Cretan(v;2; (v - k)x 2 + ky 2 ; y,x;determinant). 

We have used the notation Cretan(order;r;u.r, y,x;determinant,). □ 

In all these Cretan - Hadamard cases (but not in all cases) the Balonin-Sergeev- 
Cretan(4t- 1,2) matrix with higher determinant comes from SBIBD(4t - 1,2 t,t) while 
the SBIBD(At - 1,2t - 1, t - 1) gives a Cretan(4t - 1,2) matrix with smaller determinant. 
These examples have been given as they may give circulant SBIBD when other matrices 
do not necessarily do so. 


4. Main Equivalence Theorem 


Lemma 1 (Hadamard to SBIBD Lemma). Whenever there exists an Hadamard 
matrix of order At there exists an SBIBD{At - 1,2t - 1, t - 1). 


Proof. Let H be the Hadamard matrix of order At. We use the Hadamard equivalence 
operations to make the normalized Hadamard matrix G with the first row and column all 
+l’s. 

It now has the form 

'1 1 ... 1 ' 

1 



1 


where I? is a ±1 matrix of order 4t-1 containing 2 1 —l’s and 2 1- 1 l’s per row and column. 
B satisfies BB T = At I - J. We form A by replacing the -1 elements of B by zero, that 
is A = |( J + B). Then A satisfies AA T = tl + (t - 1) J and is the incidence matrix of an 
SBIBD{At - 1, 2t - 1, t - 1) as required. □ 


Lemma 2 (SBIBD to Cretan(4i-1,2) Lemma). Whenever there exists an SBIBD{At- 
l,2t- l,t- 1), there exists an Cretan-Mersenne 2-level matrix. 


Proof. We take the SBIBD and replace the zeros by x and the ones by y. Call this 
matrix C. 

We choose x - 1 and y - . We first show this satisfies the orthogonality equation. 

The inner product of any row with itself is 2 1 + (2 1 - 1) = u> which is a constant 

given t. Thus we have the radius equation. Using Hadamard equivalence operations we 
permute the columns of C until we have 


2t-l 


2t 


y,y,---,y,y, y, y 

y,y,---,y,y x,x,...,x,x 


00 ^ 00 j • • • ^ 00 * 00 ^ 00 j 00 ^ ^ 00 . 00 

y,y,...,y,y x,x,... y x,x 


t -1 


t 


t 


t 
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Then the inner product of rows i and j is 


(t - 1 )y 2 + 2txy + tx 2 = (t - 1) | —— j— j + 2 1 

Hence CC T = uil and we have the matrix orthogonality equation. Thus we have formed 
the required 2-level Cretan-Mersenne matrix. □ 

Lemma 3 (Cretan(4t - 1,2) to Hadamard Lemma). Whenever there exists an 
SBIBD(4t - 1, 2t - 1, t - 1), there exists an Hadamard matrix of order 4 1 . 

Proof. We take the Cretan(4t- 1,2) matrix and replace the entries which are not 1 by 
zero. This gives the incidence matrix A of the SBIBD(At - 1,2i - 1, t - 1). By definition 
of SBIBD the inner product of any pair of distinct columns is A = t -1 so it satisfies 
AA T = tl + (t - 1) J, J the matrix of all l’s. We note that it satisfies AA T =tl+(t- 1 )J. 
It has a total of 4i - 1 entries, 2t - 1 l’s and 2 1 zeros in each row and column. We now 
form B = 2 A - J, of order 4t - 1 and note BB T = At I - J. Bordering B with a row and 
column of l’s gives the required Hadamard matrix. □ 

Theorem 2 (Co-Existence of Hadamard Matrices and Cretan(4f - 1,2) 
-Mersenne Matrices Theorem). The existence of an Hadamard matrix of order At is 
equivalent to the existence of a Cretan(4t - 1, 2)-Mersenne matrix. 

Proof. We have shown that and Hadamard matrix of order At can be used to form an 
SBIBD{At - 1, 2t - 1, t - 1) using Lemma 1 as its core. 

Lemma 2 was used to show that whenever there exists an SBIBD(4t - 1,2f - 1, t - 1), 
there exists an Cretan-Mersenne 2-level matrix. 

This shows that an Hadamard matrix can always be used to obtain a Cretan-Mersenne 
2-level matrix. 

To prove the converse we use the Cretan(4t- 1,2) to Hadamard Lemma 3 to show 
how a 2-level Cretan(4t- 1,2) matrix can always be used to form an Hadamard matrix. 
Whenever there exists an SBIBD(At - 1, 2t - 1, t - 1), there exists an Hadamard matrix 
of order At . This shows that a Cretan-Mersenne 2-level matrix can always be used to 
obtain a Cretan-Mersenne 2-level matrix. 

Thus we have shown the existence of one gives the existence of the other. This 
completes the proof. □ 

Conjecture: Since Hadamard matrices are conjectured to exist for all orders At, t > 0 
an integer, 2-level Cretan(4t - 1,2)-Mersenne matrices are conjectured to exist for all 
orders At - 1, t > 0 an integer. 

5. Conclusions 

Cretan matrices are a very new area of study. They have many research lines open: 
what is the minimum number of variables that can be used; what are the determinants 
that can be found for Cretan(n; r) matrices; why do the congruence classes of the orders 
make such a difference to the proliferation of Cretan matrices for a given order; find the 
Cretan matrix with maximum and minimum determinant for a given order; can one be 
found with fewer levels? 

We conjecture that co & v will give unusual conditions. □ 
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